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, Abstract. Close-in extrasolar planets experience extreme tidal interactions with their host stars. This may lead to a reduction 

of the planetary orbit and a spin-up of stellar rotation. Tidal interactions have been computed for a number of extrasolar planets 

I in circular orbits within 0.06 AU, namely for OGLE-TR-56 b. We compare our range of the tidal dissipation value with two 

dissipation models from Sasselov <2003> and conclude that our choices are equivalent to these models. However, applied 
to the planet OGLE-TR-56 b, we find in contrast to Sasselov (.2003/) that this planet will spiral-in toward the host star in a 
few billion years. We show that the average and maximum value of our range of dissipation are equivalent to the linear and 
quadratic dissipation models of Sasselov 1 2003 1. Due to limitations in the observational techniques, we do not see a possibility 

' to distinguish between the two dissipation models as outlined by Sasselov 12003 J - OGLE-TR-56 b may therefore not serve as 

a test case for dissipation models. The probable existence of OGLE-TR-3 b at 0.02 AU and the discovery of OGLE-TR-1 13 b 

' at 0.023 AU and OGLE-TR-132 b at 0.03 AU may also counter Sasselovs <2003> assumption of a pile-up stopping boundary at 

; 0.04 AU. 
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1. Introduction 

. Since the discovery of the companion of 5 1 Pegasi (Mayor & Queloz I1995> . many more extrasolar planets orbiting sun-like 
' stars of spectral types F, G and K have been discovered in the mass range from 0.16 Jupiter masses (Mj) to 15 Mj. The total 
' number of known extrasolar planets is now 122 (June 2004), including thirteen planetary systems with two or more planets. Only 
the minimum planetary mass Mp - Mp sin / can be derived with the radial velocity method since the inclination angle / of the 
planetary orbital plane cannot be determined unambiguously. The elements of these planetary orbits are determined from the 
analysis of the variation of radial velocities of the central star. From the precise determination of the radial velocity variation 
period, interpreted as the extrasolar planet's orbital period, and an estimate of the stellar mass, the semi major axis can be derived 
' with the aid of Kepler's third law. Surprisingly, planets were found from very close distances to their host stars to distances as 
far as 4 AU (the outer known planet of 55 Cnc). So far, only one planet has allowed us to derive its true mass Mp (and size 
and therefore density) from its transit across the stellar disk of HD209458 in combination with radial velocity measurements 
(Charbonneau et al. 120001 Jha et al. 120001 Henry et al. l2000l l. In addition, upper limits of true masses exist for few other planets, 
but are not verified. 

The closest extrasolar planet discovered so far, OGLE-TR-56 b which is also a transiting planet, was reported by Konacki 
et al. J2003at from the OGLE survey (Udalski et al. 2002). Improved orbit and planetary parameters have been determined by 
Torres et al. ((20041. Sasselov (2003 1 speculates on the stability of the orbit and atmosphere. There are further possible transit 
candidates from the OGLE survey: OGLE-TR-113 b (Konacki et al.p004 Bouchy et al. 2004) and OGLE-TR-3 b (Dreizler et 
al.|2003 ) which (if the existence of the latter is verified) are both as close as OGLE-TR-56 b, OGLE-TR-10 b and OGLE-TR-58 
b reported by Konacki et al. f2003M and OGLE-TR-132 b reported by Bouchy et al. (2004i. The close vicinity of OGLE-TR-56 
b and OGLE-TR-3 b to their respective host star results in the orbit being perturbed by tidal forces exchanged between the star 
and the planet. We shall apply the equations described by Patzold & Rauer ( 2002 1 to these OGLE cases which are exact because 
of the circular nature of their orbits (e - 0). Other planets with circular orbits but which are farther away will be included in this 
work for comparison. We discuss the coupled effects of orbital decay and spin-up of the stellar rotation as a function of time in 
the context of the speculations given in Sasselov ( I2003> . 
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Table 1. Parameters of close-in extrasolar giant planets and their host stars with circular orbits with semi major axis < 0.06AU 
(as of June 2004) 



Parameter 






OGLE-TR 56 


OGLE-TR 3 


OGLE-TR-113 


OGLE-TR 10 


OGLE-TR 132 


Spectral type 






GO 


GO 


K 


GO 


F 


Planetary mass 


[Mj] 




1.450 


0.500 


1.35 


0.700 


1.010 


Semi major axis 


[AU] 




0.023 


0.023 


0.023 


0.042 


0.0306 


revolution period 


[days] 




1.212 


1.190 


1.43 


3.101 


1.690 


Orbit excentricity 






0.000 


0.000 


0.000 


0.000 


0.000 


Stellar rotation 


[days] 




18.399 


27.757" 


27.757" 


27.757" 


27.757" 


Stellar mass 






1.040'' 


1.000'' 


0.770'' 


1.000" 


1.340" 


Stellar radius 


[%] 




1.100" 


1.106'' 


0.822" 


1.106" 


1.489b 






average 


1.230 


0.870 


1.080" 


1.300'' 


1.150" 


Planetary radius"^ 


[Wj] 


maximum 


1.972 


1.383 












minimum 


0.861 


0.604 








Roche zone 


[Rq] 




2.706 


2.721 


2.022 


2.721 


3.663 


Stellar J2 


[10-«] 




-2.749 


-1.277 


-0.681 


-1.277 


-2.325 


I, 






0.07 


0.07 


0.07 


0.07 


0.07 








0.168 


0.168 


0.168 


0.168 


0.168 


System property factor 


[lO-^M^^^] 




2.290 


0.827 


0.577 


1.158 


6.386 



Parameter 






OGLE-TR 58 


BD-10 3166 


HD 209458 


HD 76700 


HD 49674 


Spectral type 






GO 


G4 


GO 


G6 


G5 


Planetary mass 


[Mj] 




1.600 


0.480 


0.690 


0.197 


0.120 


Semi major axis 


[AU] 




0.052 


0.046 


0.045 


0.049 


0.057 


revolution period 


[days] 




4.345 


3.487 


3.525 


3.971 


4.948 


Orbit excentricity 






0.000 


0.000 


0.000 


0.000 


0.000 


Stellar rotation 


[days] 




21.157' 


27.757" 


15.700 


27.757" 


27.757" 


Stellar mass 


[M,0] 




0.990" 


1.100" 


1.050" 


1.000" 


1.000" 


Stellar radius 






1.106" 


1.006" 


1.200" 


0.956" 


0.981" 






average 


1.600" 


0.850 


0.970 


0.632 


0.535 


Planetary radius*^ 


[Ri] 


maximum 




1.349 


1.540 


1.003 


0.850 






minimum 




0.589 


0.673 


0.438 


0.371 


Roche zone 


[%] 




2.721 


2.475 


2.952 


2.352 


2.414 


Stellar J2 


[10-n 




-1.290 


-0.874 


-4.855 


-0.825 


-0.892 


I, 






0.07 


0.07 


0.07 


0.07 


0.07 








0.168 


0.168 


0.168 


0.168 


0.168 


System property factor 


[lO-^M^i?^] 




2.661 


0.472 


1.676 


0.157 


0.109 



" If the true stellar rotation period is not known, the value was set, without loss of generality, to the sun's rotation period. 
" DerivedfromAUer et al. ( 19821. 

" If the planetary radius is not known, then upper and lower limits can be derived by setting the planetary bulk density to 250 ^ (Henry et 

al. l2000t Jha et al. 20Q0j and 3000 4 (Cameron et al. ll999i) . respectively. As an average, the density of 1000 was assumed. 
" This is the true radius as derived from transit observations ( Charbonneau et al.[2000 Konacki et al. l2003bV 



2. Tidal interaction 

Close-in extrasolar planets experience strong tidal interactions with their central star. If the planet's orbital period P is smaller 
than the star's rotation period P,, tidal friction will lead to a spin-up of the star and, due to the conservation of momentum, will 
also lead to a decrease of the semi major axis of the planet's orbit.The parameters of the planets and the stars considered in this 
work are listed in Table ^ Only planets with circular orbits within 0.06 AU are used in this study. A more detailed study for 
elliptical orbits in general is in preparation. 

The rotation periods of only two stars (OGLE-TR-56 via vsin/ and HD209458) in Tableware known. Due to the small 
revolution periods of the planets, all other stars must be fast rotators (P» < 3days) and the condition P» > P for the decrease of 
the orbit does not hold and is considered as not very probable. The stellar rotation periods are set to the solar rotation period of 
27 days if these parameters are not explicitly known or estimated otherwise. Following Murray & Dermott (1999) , Goldreich & 
Soter J1966> . Goldreich <& Nicholson <1977> and Zahn (I1977> . we derive for the change in stellar rotation and the change in orbit 
radius 

— = -sign(Q.-«)^^-J^-3j-^ (1) 
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Fig. 1. a) The time scales for spiralling into the host stars (according to equation (|5}) of the extrasolar planets located at semi 
major axis of today. The symbols mark the avarage value of j^, the error bar the selected range of ^ (see section|2li. b) The 
Doodson constant of extrasolar planets with e - Q for a semi major axis less than 0.06 AU. 



da . 3k2jMp\(R,\^ IGM, 
— = sign(i2«-n) — — fl-i/ — — 



(2) 



where Q, and n are the stellar rotation and planetary revolution, respectively, k2t is the stellar Love number, C, is the stellar 
moment of inertia, g* is the stellar dissipation factor which includes tidal friction, M» and are the stellar mass and radius, 
respectively, a is the semi major axis of the planetary orbit and G is the gravitational constant. The relations Q and (|3 are 
valid for a circular orbit (e = 0) which is assumed from the observations for the OGLE candidates. We used the stellar moment 
of inertia C* which is usually expressed in Q as C* = /«M*7?J where /* is the normalized stellar moment of inertia. The mass 
distribution of the standard solar model, as well as the GONG model, yield an /» = 0.07 (Bahcall et al. 120011 Christensen- 
Dalsgaard et al.'1996'l. Both the stellar Love number ^2* and the tidal energy dissipation factor Q, are only poorly known. We use 
10*' < < 3 ■ 10^ and 0.02 < /t2. < 0. 17 which results in 5 ■ 10*' < < 1 .5 ■ 10'' (Patzold & Rauer 2Q02 ) for these two combined 

parameters, with ^ — 1.2 ■ 10^ as an average. This is in agreement with the range of energy dissipation for main sequence stars 
used in the literature (Lin et al. 1996 Trilling et al. 119981 Bodenheimer et al. 2001 1. Fixing the Love number at ^2* = 0.17 
(Patzold & Rauerl!2QQ2t, the average dissipation factor is in agreement with the value found by Patzold & Rauer (2002^ for the 
F-stars; a Love number of k2, - 0.02 yields - 2.4 ■ 10^. 

The stellar radius is another important input parameter for the relations Q and (|2j. We derive these values from calibration 
curves based on stellar evolution models (Aller et al. ll982> . Integration of Q yields for the stellar rotation 



a(f) = - 



af + sign(Q, - n)^4^^Rl 



2 e. 



1/13 



(3) 



and for the semi major axis 



a(t) 



+sign(Q-«)y-^^/?5^/^ 



(4) 



The time scale for spiralling inward to the stellar Roche limit Oioche - 2.46 (Chandrasekhar [T^87> from the present semi major 
axis fl() is (as derived from @) 



(5) 



Here and in all following equations we have replaced the true planetary mass Mp with the minimum mass Mp - Mp sin /. Aside 
from the larger uncertainties of Q,, the time scale To, in Q has to be considered as an upper limit because the true planetary 
mass (Mp > Mp) will further decrease this value. The time scale Ta* represents the remaining time for the tidal process and not 
the time since formation of the stellar system. 



4 



M. Patzold et al.: Tidal interactions of close-in extrasolar planets: the OGLE cases 



O/kJ' = l.DeS average ^ 
QM- 1 5.0e6 minnnjm Y 




0.01 



0.03 



3 
< 



K 



i? 
E 

E 
» 

10 



0.OT 




0.0 0.5 1.0 1.5 2.0 2.S 3.0 3.5 4.0 4.S S.O 
time [10^ years] 



0.0 0.S to 1.5 2.0 25 3.0 3.S 4.0 4.S S.O 
time [10^ years] 



(b) 



Fig. 2. The change in semi major axis of OGLE-TR-56 b a) and OGLE-TR-3 b b) within 5 bilHon years for different values 
For values of ^ smaller than the avarage, the planet needs less than a few billion years to spiral into the star 



3. Discussion 

Figure [2 (a) shows the time needed to let the planets spiral into the central star starting from their orbital position of today. The 
error bars reflect the variation of the parameter The closest planets, OGLE-TR-56 b and OGLE-TR-3 b, will spiral in within 

a few billion years on average, even much faster considering the lowest boundary of The upper boundary of ^ is obviously 
generously chosen so that no influence on the orbit can be seen. This orbit may be considered as stable. Possible candidates for 
strong orbital decay can be found by examining the system property factor -^^Rl (Patzold & Rauer 2002 1. Another possibility 
is the well known Doodson constant (Figure [2 (b)) which describes the amplitude of tidal forces exchanged between the two 
partners: 



(6) 



Although our model does include viscosity and friction inside the star expressed by the dissipation factor 2*> we do not agree 
with Sasselov ( 2003 ) who considers the existence of OGLE-TR-56 b as counterproof of the orbit decay description in Patzold & 
Rauer (2002). The computed time scale Q is counted from today to the time of planetary disruption and not from the time of 
creation of the planetary system. Having in mind the many theories of planet formation, migration and tidal movement, it is very 
difficult to assess how the close-in planets finally achieved the current observed positions. Let us compare the equation Q used 
by Patzold & Rauer ilQOli and the relation used by Sasselov (;2003t . which was given by Rasio et al. Gl996i) : 



fl 1 / Mcz Mp 



1 + 



Mp 



a 



(7) 



where / - min (;;^)] or / = min |^1, (5^)] is the total integrated turbulent viscosity in the quadratic viscosity suppression case 
and the linear suppression case, respectively, tc - 18 days is the turbulent time scale (Sasselov 2003 1 and Mcz - O.IM, is the 
mass of the convection zone for a sun-like star as derived from the standard solar model. Comparing directly ^ with Q, we 
have: 

= I ^ for the quadratic suppression case and 
^ = l^'^l^c^ for the linear suppression case. 



The former one yields ^ 2.7 ■ 10'" which is comparable to the upper boundary of our chosen range of while 

-* i quadratic 

the latter one yields ^9-10^ which is slightly smaller than our chosen average value of Q,. Figure |2 shows the time 

-* I linear 

dependence of the orbital decay of OGLE-TR-56 b and OGLE-TR-3 b. The three curves describe the minimal, average and 
maximal value of ^ as discussed above in (|2j. For OGLE-TR-56 b, the average value of ^ = 10*^ lets the planet spiral into the 

star (Roche limit of 2.46/?*) within a few bilhon years, the lower value of £^ achieves the same result in about 1 .5 bilhon years. 
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Fig. 3. The change in stellar rotation period of the host stars of OGLE-TR-56 b (a) and OGLE-TR-3 b (b) due to tidal friction. In 
this case, it was assumed that the stellar body is rigid. For values of ^ smaller than the avarage, the star is significantly spun up. 

If the linear suppression case which considers dissipation is true, the planet is doomed, since in that case ^ is even smaller than 
our chosen average value. No effect within the life time of the star is observed for the upper boundary value of ^ which is more 
or less identical with the quadratic suppression case. Figure |2(b) shows the temporal dependence of OGLE-TR-3 b. Although 
this planet is slightly further away from its central star than OGLE-TR-56 b, its mass is lower (0.5Mj (Dreizler, priv. Comm.)) 
and its Doodson constant is smaller The result is that the tidal effects are weaker than for OGLE-TR-56 b (Figure [2(b)) and 
FigureEKb)). 

Due to the conservation of angular momentum, the star is spun up as the orbit decays. Equation Q describes the stellar 
spin-up by assuming the star being a solid body, represented by the introduction of the total moment of inertia relative to the 
rotation axis C*. We neglected friction within the star. Figure |3] shows that the spin-up of the entire solid body according to ([0 
leads to a rapid decrease of the rotation period from 18 days to less than 15 days within the time of 800 Million years needed for 
the planet to spiral-in towards the stellar Roche zone. For the upper boundary of ^ no significant effect can be seen. It may be 
argued that mass loss of the star due to stellar wind decreases the angular momentum with increasing stellar age and thus leads 
to a decrease in stellar rotation, counteracting the tidal spin-up. But we assume that this effect can safely be neglected since our 
dataset contains only sun-like stars which will not lose a considerable amount of mass during the next four to five billion years. 
However, we did not consider friction within the stellar body and effects of the stellar magnetic field which might counteract the 
tidal spin-up of the stellar body. 

There is discussion that not the entire stellar body is spun up, but only the convection zone (Rasio et al. 119961 Sasselov 
12003 V That would mean that in Q the moment of inertia of the entire body needs to be replaced by the moment of inertia of the 
convection cell. 



— — = -sign (12, - n) 



(M, 



dt " 2CczG. 



Rl\GM (8) 



,3 



where Ccz = | Zi,=,(, '""i - ^czM^RI and Iqz is the normalized moment of inertia of the convection zone, is the 

radius of the inner boundary of the convection zone in stellar radii, and m, is the mass of the layer between r,+i and r, in stellar 
masses. Using again the standard solar model, the normalized moment of inertia of the convection zone is Iqz = 0.01, nicz - 0.01 
and Tc = 0.75. Figure|4]shows the change of the rotation of the stellar convection zone. 

4. Conclusion 

We do not agree with Sasselov <2003> that his derived long orbit decay time of about 14.3 billion years for the quadratic sup- 
pression case is a proof of his description of turbulent viscosity dissipation, simply justified by the fact that the planet still exists. 
The derived time scales in Q describe the remaining time for orbit decay and not the time since the formation of the system. On 
average, OGLE-TR-56 b will spiral into the Roche zone of its respective host star within 3.5 billion years. Because of its lower 



mass, OGLE-TR-3 b may be safe, if this planet exists at all. Correcting Sasselov ( I2003> . we have shown that the suppression 
case for dissipation in the stellar convective layers is considered by our choice of Assuming the linear suppression case, the 
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Fig. 4. The change in stellar rotation period of the host stars of OGLE-TR-56 b a) and OGLE-TR-3 b b) due to tidal friction. In 
this case, it was assumed that the tidal enery is dissipated in the outer stellar layers. Neglecting friction between the outer layers 
and the core, the spin-up of the outer stellar layer is significant even for the maximum value of 



dissipation in the convection zone is considered as our average value of j^. The quadratic suppression case is covered by our 
upper boundary value for Applying this model would yield absolutely no tidal effect for any planet at any distance. 

Would it be possible to distinguish between these two dissipation models? This would mean to measure the orbital decay 
rate of the planet via the transit time across the stellar disk. The velocity in a circular orbit is = = Assuming that the 
orbital radius a is decaying with the rate at first order a{t2) = a{ti) - ^{h -h) it is clear that 



, p/M 



Ri da, , 3k2. Mp (R, 

' ' ) ^g¥^''-''^ 

Q. 



Thus, the transit times would change by 0. 1 seconds in one hundred years or by 1 ms per year for the average value of , 
and would require a timing accuracy of the transit times better than 0.1 millisecond. With the current observation methods 
and integration times this accuracy is simply not achievable. Therefore, one would not be able to judge which of these models 



describes the convective dissipation and the proposal of Sasselov J2003> to use OGLE-TR-56 b as a test case for convection 
dissipation models cannot hold. 

We have also shown that the consequence of orbital decay results in the spin-up of the stellar rotation. For all selected values 
of Y", the convection zone is spun-up significantly. If this is true, one should find all stars with close-in extrasolar planets to be 
fast rotators. However, this is not observed. Friction between the convection zone and the remaining inner radiation zone and 
core avoids a fast acceleration of the convection zone, in particular of the observable outer layers. 

We also do not agree with the conception of a 'pile-up' of close-in planets at about 0.04 AU due to an unknown stopping 
mechanism for the inward migration (Sasselov l2003t . First, this concept is based on an extremely limited observation base which 
does not allow any conclusions. Second, we have shown here that the orbital decay of planets closer than 0.04 AU is fast if the 
planetary mass is sufficiently large and only in rare cases allows the observation of planets within 0.04 AU which are in the 
process of spiraling toward the host star. Third, OGLE-TR-3 b, if it exists at all, is the counterproof: due to its low mass the tidal 
interaction is negligible. 
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